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I. INTRODUCTION 

Hadrons are the smallest material objects in the uni- 
verse known to be of finite size. The building blocks of the 
Standard Model, including leptons and quarks, are not 
known to have finite size with the smallest scale set by 
the experimental limit on the size of the electron, which 
is smaller than lCU 18 m in diameter [TJ. Hadrons are dis- 
tinguished in two families: mesons, which are made out 
of a quark and an antiquark pair and baryons, which are 
made out of three quarks. Quantum Chromodynamics 
(QCD) does not exclude the possibility of other forms of 
hadronic matter such as di-baryons or pentaquarks, but 
none of these have been found thus far. The typical scale 
of a hadron radius is set by the well known charge radius 
of the proton, which is equal to 0.8768 (69) x 10~ 15 m [2]. 
The very concept of size, both classically and quantum 
mechanically, raises the issue of shape and it is therefore 
natural to inquire about the shape of hadrons. The shape 
of hadrons concerns microscopic objects at the scale of a 
femtometer (10~ 15 m). 

Inquiring about the shape of a subatomic particle is 
equivalent to raising the question whether the distribu- 
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tion of its constituents or some of the extensive and there- 
fore distributed properties, such as mass or charge, de- 
viate from a spherical distribution, which is assumed to 
be the default distribution. However, hadrons are objects 
which are understood only within a quantum mechanical, 
relativistic framework. It is thus necessary to reexamine 
critically the concept of size and shape in the context of 
quantum mechanics and relativity before we address the 
question of the shape of hadrons. In parallel, we will also 
address the issue of how sizes and shapes are determined 
for particles of the microcosm. 

Knowledge of the shape of the fundamental building 
blocks of the Universe, is not a curiosity, although it cer- 
tainly comes close to being an example of the Aristotelian 
claim of the intrinsic human need to "know" . Experience 
from the determination and subsequent understanding of 
shapes of other objects in the microcosm such as those 
of atoms and nuclei, shows that this line of investigation 
is particularly fertile for the understanding of the inter- 
actions of their constituents amongst themselves and the 
surrounding medium. For hadrons this means the in- 
terquark interaction and the quark - gluon dynamics. 

While the theoretical foundations describing hadrons 
as the smallest objects in the universe to which size and 
shape can be attributed are solid, the empirical knowl- 
edge concerning shape is limited and derives only from 
the detailed study of the transition to the first excited 
state of the proton, the A + (1232) resonance [3 . It is 
interesting to observe that while the determination that 
hadrons have size emerged early on, through the seminal 
work of R. Hofstadter and collaborators [4] and played 
a leading role in guiding hadron research ever since, the 
determination of shape has been elusive and continues to 
be very limited. 

In the rest of this section we will review the develop- 
ment of the concepts of size and shape in classical and 
quantum mechanics with and without relativity, so as to 
establish the appropriate language needed to discuss the 
topic of "Shape of Hadrons" . This is necessary as hadrons 
are systems requiring a relativistic quantum mechanical 
description. In doing so, we will provide some historical 
background on how these concepts have developed. 



A. Historical Development 

The issue of shape of subatomic particles arose most 
acutely in the case of nuclear physics. It is interest- 
ing to observe that historically the issue concerning the 
shape of atoms being non spherical never caused much 
surprise, perhaps because of the planetary (Rutherford) 
model, which intrinsically invokes non-spherical shapes. 
The discovery by Rabi and collaborators [5] that the 
deuteron had a static quadrupole moment and therefore 
its shape was not spherical was regarded as a major sur- 
prise. The discovery of deformation in the deuteron, and 
nuclei in general, was interpreted correctly as arising due 
to the existence of non-central (tensor) forces among nu- 



cleons. Shortly afterwards, Gerjuoy and Schwinger pro- 
posed that trinucleon deformation (e.g. 3 H) could re- 
solve some peculiarities in the spectroscopy of those sys- 
tems [5J. This conjecture proved to be wrong - the ef- 
fects were eventually understood to be due to mesonic 
degrees of freedom. The deuteron and trinucleon cases 
dramatically showed that understanding the shape of a 
subatomic particle requires a detailed knowledge of its 
constituents and it provides important information for 
their dynamics. Following the success of Rabi, the es- 
tablishment and quantification of deformation through 
the measurement of the electric quadrupole moment of a 
particle was widely employed to map the systematics of 
deformation in atomic nuclei. 

Unfortunately a number of misconceptions arose as 
a result of the successful use of the determination of 
quadrupole moments in inferring deviations from spher- 
ical shapes for atomic nuclei. The fact that the mea- 
surement of a quadrupole moment is possible only for 
systems (particles) possessing spin equal to one or bigger 
led incorrectly to the belief that the shape of particles 
possessing spin or 1/2 cannot be determined. The im- 
possibility to measure a quadrupole moment of such par- 
ticles was mistakenly interpreted as signifying a spherical 
shape. It took more than two decades before this issue 
was clarified, primarily through the work of P. Brix and 
coworkers [7]. The realization that the determination of 
the intrinsic shape of a system is quite distinct from the 
ability to measure its quadrupole moment helped the field 
develop. Nuclear physicists developed new techniques to 
measure shapes that were also applicable to nuclei of spin 
or 1/2. In cases where the rigid rotor (shape) approx- 
imation could be made, their excitation spectrum could 
be used to reconstruct the density [8]. 




r(fm) 

FIG. 1: Tomographic view of the deformed nucleus 154 Gd 
derived from the study of the rotational bands of this nucleus 
using electron scattering [8]. The plot shows the contours 
of equal charge density revealing its shape. This spin zero 
nucleus has a vanishing quadrupole moment and therefore 
what is revealed through this reconstruction is its intrinsic 
shape. 
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FIG. 2: The transition densities characterizing the moving 
charges involved in excitations in the deformed nucleus 154 Gd, 
derived from the study of this nucleus using electron scatter- 
ing [S] , demonstrate that vibrations along the two axes of the 
ellipsoidal shape have different spatial extent. 



In the seventies and eighties high resolution electron 
scattering, using this technique was used to map the 
shapes of the deformed nuclei in the rare earth and ac- 
tinidc regions of the periodic table. A superb example 
of the finesse of this tool is shown in Figs. [T]and[2] The 
reconstructed ground state charge density of the 54 Gd 
spin zero nucleus [8] is mapped with high precision and 
the isodensity contours reveal vividly its deformation, as 
is seen in Fig. [I] Similarly the transition densities of 
154 Gd, shown in Fig. [2] provide a vivid geometric repre- 
sentation of the shape oscillations of the system along 
the long and short axes of symmetry again revealing its 
non-spherical charge distribution [5]. 

The understanding of the shape of nuclei led to a better 
understanding of nuclear structure and an appropriate 
language to describe a number of important nuclear phe- 
nomena. The development of the formalism of shape os- 
cillations viewed as normal modes of the oscillating quan- 
tum liquid, nuclear matter, was a crucial milestone in the 
field of nuclear physics and indeed of physics. Modes of 
shape oscillations such as the "giant dipole oscillation", 
the "breathing mode" or the "scissor's mode" [5j have 
yielded valuable knowledge on a number of parameters 
characterizing nuclear matter. For instance, its com- 
pressibility a parameter of critical importance in the un- 
derstanding of supernova explosions, is derived from the 
study of the breathing mode of nuclei. Nuclear shapes 
and shape oscillations have also led to paradigms, which 
are driving the development in other fields of physics, 
such as the observation of a "scissor's mode" in Bose- 
Einstein condensates in low temperature physics |10l lllj . 

It may appear from the preceding introductory com- 
ments that to inquire about the shape of an object pos- 
sessing size is an obvious undertaking. However, it took 
more than twenty five years from the indication of the fi- 
nite size of the proton to the inquiry about its shape. The 



conjecture that hadrons would have non-spherical ampli- 
tudes was first made by Glashow in 1979 on the basis 
of non-central (tensor) interactions between quarks |12j . 
Glashow argued that this would resolve a number of in- 
consistencies that QCD was facing at the time if the con- 
strain of sphericity of the shape of hadrons was relaxed. 
The conjecture of non-spherical hadrons originally was 
based on the premise that there is a color spin-spin in- 
teraction between the quarks [13] , which is modeled after 
the interaction between magnetic dipoles in electromag- 
netism, the so-called "Fermi-Breit" interaction [TJ. A 
few years later Isgur, Karl, and Koniuk wrote a semi- 
nal paper j!5) . which offered an impressive list of indi- 
rect empirical evidence for this hypothesis. However, in 
a remarkable similarity to the flawed trinucleon defor- 
mation hypothesis of Schwinger, due to the oversimpli- 
fied description of the system, the non-relativistic shell 
model description of baryons ("tri quarks") is now also 
found to be unable to quantitatively describe the defor- 
mation when solely invoking the color magnetic tensor 
interaction. The inadequacy of the non-relativistic de- 
scription and of the phenomenological description of the 
constituents used (lack of mesonic degrees of freedom) are 
understood to be the principal deficiencies of this model. 

In their paper, Isgur, Karl, and Koniuk singled out the 
quadrupole amplitude in the A — > Nj transition as be- 
ing a most sensitive test of this hypothesis. Of additional 
interest arc the quark model calculations, which showed 
that the D-state admixtures caused by the color hypcrfinc 
interaction predict a non-zero neutron charge distribu- 
tion and root-mean-square (RMS) charge radius [T5HIH]- 
These theoretical speculations induced concerted exper- 
imental and theoretical efforts to measure and calculate 
deviations from spherical symmetry (non-spherical am- 
plitudes) in hadrons. 



B. Size and Shape in Classical and 
Non-Relativistic Quantum Mechanics 

The concepts of both size and shape, because of their 
familiarity in everyday language, are often taken to be 
intuitively apparent, at least in classical physics. How- 
ever, a careful examination reveals that this is not at all 
the case except for rigid objects with uniform density and 
sharp boundaries. The size of a hurricane or the size and 
the shape of nebula (e.g. the crab nebula) are not easy to 
quantify. However, the distribution in space and time of 
some extensive property of an object such as its mass or 
charge can uniquely and unambiguously be defined. Its 
mass density p(r) is uniquely defined and so is its vari- 
ation in time p(r,i). In classical physics densities can 
be precisely defined and measured and their knowledge 
allows one to define a "size" and a "shape". Moments 
of the density distribution are often quoted, which, in 
simple geometrical limiting cases, have the expected cor- 
respondence to the naive concept of size or shape. For 
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instance the second moment of the density distribution 

(r 2 ) = J dr ■ r 2 ■ p(r) (1) 

corresponds to the radius of a spherical body with uni- 
form distribution p(r) = po9(R — r). 

For objects whose density distribution deviates from 
spherical symmetry it is obvious that higher moments 
will assume non- vanishing values. The first such moment 
whose non vanishing value indicates non-sphericity is the 
quadrupolc moment Qij : 

Qij = (Qij) = Jdr ■ (3r irj - r%j) • p(r), (2) 

with i,j = 1,2,3 denoting the spatial directions. It is 
worth noting that it is possible to have non-spherical dis- 
tributions that have vanishing quadrupole moments. 

The introduction of non-relativistic quantum mechan- 
ics and the implications of the uncertainty principle have 
influenced profoundly our understanding of the concept 
of size. Early on, with the aid of the "correspondence 
principle" it was realized that "size" could be expressed 
in terms of the RMS radius given by Eq. ([T]) where 
p(r) = (r)?/>(r) is the probability density expressed in 
terms of the wave function ^(r) of the object. Likewise, 
the quadrupole moment of a system (Qij) manifests devi- 
ation of its probability density from spherical symmetry. 

C. Size and shape in relativistic systems 

The introduction of relativity does complicate matters. 
It is well understood that both the size and the shape of 
an object, are not relativistically invariant quantities: ob- 
servers in different frames will infer different magnitudes 
for these quantities. Furthermore when special relativity 
is written in a covariant formulation, the density appears 
as the time (zeroth) component of a four-current density 
= (p, J) (in units where the speed of light c = 1). 

Besides the relativistic kinematical effects, e.g. due 
to length contraction, the concept of size and shape in 
relativistic quantum systems, such as hadrons, is also 
profoundly modified as the number of constituents is not 
constant as a result of virtual pair production. Consider, 
as an example, a hadron such as the proton, which is 
probed by a space-like virtual photon, as shown in Fig. [3] 
A relativistic bound state is made up of almost massless 
quarks. The three valence quarks, which make up for the 
proton quantum numbers, constitute only a few percent 
of the total proton mass. In such a system, the wave 
function contains, besides the three valence quark Fock 
component \qqq), also components where additional qq 
pairs, so-called sea-quarks, or (transverse) gluons g are 
excited, leading to an infinite tower of \qqqqq), \qqqg), ■■■ 
components. When probing such a system using electron 
scattering, the exchanged virtual photon will couple to 
any quark, both valence and sea in the proton as shown 




FIG. 3: Coupling of a space-like virtual photon to a relativis- 
tic many-body system, as a proton. Upper panel : diagonal 
transition where the photon couples to a quark, in the leading 
3q Fock component (left), or in a higher 5q Fock component 
(right). Lower panel : process where the photon creates a qq 
pair leading to a non-diagonal transition between an initial 
3q state and a final 5q state in the proton. 



in Fig. [3] (upper panel). In addition, the virtual photon, 
can also produce a qq pair, giving rise e.g. to a transition 
from a 3q state in the initial wave function to a 5q state in 
the final wave function, as shown in Fig. [3] (lower panel) . 
Such processes, leading to non-diagonal overlaps between 
initial and final wave functions, are not positive definite, 
and do not allow for a simple probability interpretation 
of the density p anymore. Only the processes shown in 
the upper panel of Fig. [3] with the same initial and final 
wave function yield a positive definite particle density, 
allowing for a probability interpretation. 

This relativistic dynamical effect of pair creation or an- 
nihilation fundamentally hampers the interpretation of 
density and any discussion of size and shape of a rela- 
tivistic quantum system. An interpretation in terms of 
the concept of a density requires suppressing the contri- 
butions shown in the lower panel of Fig. [3j This is pos- 
sible when viewing the hadron from a light front, which 
allows to describe the hadron state by an infinite tower 
of light-front wave functions. Consider the electromag- 
netic (e.m.) transition from an initial hadron (with four- 
momentum p) to a final hadron (with four-momentum 
p') when viewed from a light front moving towards the 
hadron. Equivalently, this corresponds with a frame 
where the hadrons have a large momentum component 
along the z-axis chosen along the direction of the hadrons 
average momentum P = (p+p / )/2. One then defines the 
light-front plus (+) component by a = a ± a 3 , which 
is always a positive quantity for the quark or anti-quark 
four-momenta in the hadron. When we now view the 
hadron in a so-called Drell-Yan frame |19j , where the vir- 
tual photon four-momentum q satisfies q + = 0, energy- 
momentum conservation will forbid processes where this 
virtual photon splits into a qq pair. Such a choice is pos- 
sible for a space-like virtual photon, and its virtuality is 
then given by q 2 = —q± = —Q 2 < 0, where <fj_is the 
transverse photon momentum (lying in the xy-pane). In 
such a frame, the virtual photon only couples to forward 
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moving partons, i.e. only processes such as those shown 
in the upper panel in Fig.|3]are allowed. We can then de- 
fine a proper density operator through the + component 
of the four-current by J + = J° + J 3 [20 . For quarks it 
is given by 



J+ 



qj + q 



2q\q+, with g + = (1/4)7-7+3, (3) 



where the q+ fields are related with the quark fields q 
through a field redefinition, involving the ± components 
of the Dirac 7-matrices. The relativistic density opera- 
tor J + , defined in Eq. ([3| , is a positive definite quantity. 
For systems consisting of u and d quarks, multiplying 
this current with the quark charges yields a quark charge 
density operator given by J+(Q) = + § u(0)j + u(0) — 
| J(0)7 + (0)d(0). Using such quark charge density oper- 
ator, one can then define quark (transverse) charge den- 
sities in a hadron as 21, 22] : 



Px(b) 



d 2 q± 



-iq±-b 



1 



2P^ 



(4) 



x(P 4 



qx 
2 



<i±. 



A|J+(0)|P+,-^,A), 



with A the hadron (light-front) helicity. In the two- 
dimensional Fourier transform of Eq. Q, the two- 
dimensional vector b denotes the quark position (in 
the ccy-plane) from the transverse center-of-momcntum 
(cm.) of the hadron. It is the position variable conju- 
gate to the hadron relative transverse momentum q±. 
The quantity p\(b) has the interpretation of the two- 
dimensional (transverse) charge density at distance b = 
\b\ from the transverse cm. of the hadron with helic- 
ity A. In the light-front frame, it corresponds with the 
projection of the charge density along the line-of-sight. 




J + in eigenstates of transverse spin 



as 



PTs ± ( b ) 



(27T) 2 



-iq±-b 



1 



2P 



(5) 



x(P+, q ±, s± \J+\P+, 



where s± is the hadron spin projection along the direction 
of S± . Whereas the density p\ for a hadron in a state of 
definite helicity is circularly symmetric for all spins, the 
density pr Si depends also on the orientation of the posi- 
tion vector 6, relative to the transverse spin vector Sj_, as 
illustrated in Fig. |4j Therefore, it contains information 
on the hadron shape, projected on a plane perpendicular 
to the line-of-sight. The matrix elements of the density 
operator can be written in terms of helicity amplitudes 
which in turn can be expressed in terms of the form fac- 
tors. From Ptsj_j one can then straightforwardly define 
e.m. moments quantifying the shape. As an example, 
for a hadron with spin 5 > 1/2, and with transverse spin 
orientation 5j_ = e x , the electric quadrupole moment is 
given by : 



d 2 b(bl~b 2 v )p Ts± (b) 



(6) 



These light-front densities require us to develop some 
new intuition, as they are defined at equal light-front 
time (x + = 0) of their constituents. When constituents 
move non-relativistically, it does not make a difference 
whether they are observed at equal time (t = 0) or equal 
light-front time (x + — 0), since the constituents can only 
move a negligible small distance during the small time 
interval that a light-ray needs to connect them. This is 
not the case, however, for bound systems of relativistic 
constituents such as hadrons [251 For the latter, the 
transverse density at equal light-front time can be inter- 
preted as a 2-dimensional photograph of a 3-dimensional 
object, reflecting the position of charged constituents at 
different times, which can be (causally) connected by a 
light ray. 



FIG. 4: Schematic view of the projection of the charge den- 
sity along the line-of-sight (perpendicular to the figure), for a 
hadron polarized along the direction of S±. The position of 
the (quark) charge inside the hadron is denoted by b. 

The quark charge densities in Eq. Q do not fully de- 
scribe the e.m. structure of the hadron, e.g. for spin 1/2 
the densities with A = ±1/2 yield the same information. 
We do know however that there are two independent e.m. 
FFs describing the structure of the nucleon. In general, a 
particle of spin 5 is described by (25+ 1) e.m. moments. 
To fully describe the structure of a hadron one also needs 
to consider the charge densities in a transversely polar- 
ized hadron state, denoting the transverse polarization 
direction by 5j_. The transverse charge densities can be 
defined through matrix elements of the density operator 



II. MEASURING AND CALCULATING THE 
SHAPE OF HADRONS 

The determination of the shape of hadrons, interesting 
as it may be, presents a particularly difficult situation 
both theoretically and experimentally. The challenge lies 
in identifying the observables that can provide a charac- 
teristic signal, which can be experimentally accessed with 
sufficient accuracy and can be interpreted reliably to ex- 
tract the information about shape. This has proved to be 
a particularly hard task for a number of reasons, which 
are discussed in this section. 

It has been possible in the last decade to reach the ap- 
propriate sensitivity and technical maturity to obtain and 
analyze the data that can provide the first convincing in- 
formation on the shape of hadrons. To interpret the data 
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in terms of hadronic structure quantities requires a reli- 
able reaction framework. Such a reaction framework, as 
well as the interpretation and its connection to QCD, pri- 
marily through lattice gauge calculations, have advanced 
to maturity in recent years. 

In this section we review and present these advances, 
the experimental methods, and the theoretical frame- 
work, which have allowed the first determination on the 
shape of hadrons. 



A. Empirical information for spin-1 particles: W 
boson and deuteron 

We start by discussing the empirical information on the 
c.m. moments of spin-1 particles, which are the particles 
with the smallest spin where a quadrupole moment can 
be measured. In nature, charged spin-1 particles include 
the W gauge bosons in the Standard Model of particle 
physics, the vector mesons in hadronic physics and the 
deuteron in nuclear physics. For a spin-1 system, it is 
customary to denote the three elastic e.m. form factors 
(FFs) as measured in elastic electron scattering by Gc 
(Coulomb monopole), Gm (magnetic dipole), and Gq 
(Coulomb quadrupole), where the multipole nomencla- 
ture refers to a Breit frame interpretation. 

From the empirical knowledge of the spin-1 FFs, one 
can map out the charge densities in a spin-1 particle 
of transverse polarization by working out the Fourier 
transform in Eq. ([5]), which yields monopole, dipole and 
quadrupole field patterns in the charge density [23] . The 
monopole field pattern corresponds to a circularly sym- 
metric two-dimensional distribution for a spin-1 particle 
of fixed helicity. The dipole field pattern in the charge 
distribution is specific for a relativistic theory. Indeed, a 
magnetic dipole moment in a rest frame manifests itself 
as an electric dipole moment when seen by a moving ob- 
server, proportional to the vector product (velocity) x 
(magnetic moment). The induced electric dipole moment 
(EDM) corresponding to the transverse charge densities 
Pt s± of Eq. ([5]) for transverse spin projections sj_ = 0, ±1 
is given by : 

d s± = e J d 2 bbp Ts± (b). (7) 

For example, when the transverse spin projection s± = 1, 
the expression for the electric dipole moment is [23] 

di = - (£l x e 2 ) [G M (0) - 2] e/(2M), (8) 

where M is the mass of the particle. Expressing the 
spin-1 magnetic moment in terms of the ^-factor, i.e. 
Gm(0) = g, one sees that the induced EDM d\ is pro- 
portional to g — 2. The same result was found for the 
case of a spin- 1/2 particle [23]. One thus observes that 
for a particle without internal structure, corresponding 
with g = 2 at tree level [2S], there is no induced EDM. 



The electric quadrupole field pattern in the transverse 
charge density Pts± yields a quadrupole moment, which 
is obtained, for s± = 1, from Eq. ^ as [23] : 

Qi = (1/2) [(G M (0) - 2) + (G Q (0) + 1)] (e/M 2 ). (9) 

For a charged spin-1 particle without internal structure, 
exemplified by the W gauge bosons of the standard elec- 
troweak theory, it is required that at tree level Gm(0) = 2 
and Gq(0) = —1. For elementary particles, any de- 
viations at tree level from these values would indicate 
new, beyond standard model, physics, and will show up 
in the presence of anomalous WW-y couplings, usually 
parametrized in terms of two new couplings k 7 and A 7 , 
appearing in an effective Lagrangian. In terms of those 
parameters, the W magnetic dipole and quadrupole mo- 
ments take on the values [2S] : 

fi w = e/(2M w ) [2 + (k 7 - 1) + A 7 ] (10) 
Q w = -e/M^ [1 + (/t 7 - 1) - A 7 ] (11) 

with M\y the PF-boson mass. The Standard Model val- 
ues Gm(0) = 2 and Gq(0) = — 1 equivalently corre- 
spond with k 7 = 1, A 7 = at tree level. The mea- 
surement of the gauge boson couplings and the search 
for possible anomalous contributions due to the effects of 
new, beyond Standard Model, physics have been among 
the principal physics aims at LEP-II. They have been 
prominently studied in the e + e~ — > W + W~ process 
through an s-channel virtual photon exchange mecha- 
nism. The most recent Particle Data Group (PDG) fit 
for the anomalous WWj couplings based on an analysis 
of all LEP data is given by [3D] : 

K7 = 0.973iS;gJ* A 7 = -0.0281°;"?. (12) 

One thus sees that present day information shows no ev- 
idence for anomalous WWj couplings, confirming the 
point particle values Gm(0) = 2 and Gq(0) = — 1 for the 
W bosons, leading to vanishing induced electric dipole 
and quadrupole moments according to Eqs. ([8j [9]). It is 
thus interesting to observe from Eq. ^ that Q s± is only 
sensitive to the anomalous parts of the spin-1 magnetic 
dipole and electric quadrupole moments, and vanishes for 
a particle without internal structure. 

For composite particles, it is the deviation from these 
benchmark values that indicate deformations of the 
states. A well studied example of a nuclear state is 
the deuteron. Its magnetic dipole moment is given by 
Gj f (0) =1.71 [2], close to a spin-1 particle's natural (i.e 
point-like) value. However, in contrast to the W gauge 
bosons, the deuteron has a large anomalous quadrupole 
moment. Its measured value is G^(0) = 25.84±0.03 [31] . 
Its large value was interpreted to arise from the promi- 
nent role of the one-pion exchange tensor interaction. 
One also sees from Eq. ^ that the natural value Gq (0) = 
— 1, arising in a relativistic quantum field theory for a 
spin-1 point particle, only amounts to a few percent of the 
deuteron's total quadrupole moment. For an understand- 
ing of its static properties, the deuteron can therefore be 
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FIG. 5: Two-dimensional charge density pri, according to 
Eq. for a deuteron polarized along the positive x-axis. 
The light (dark) regions correspond with largest (smallest) 
densities. The density is calculated from empirical informa- 
tion for the deuteron e.m. FFs [32]. Figure from Ref. [24] , 

considered, to a good approximation, as a non-relativistic 
bound state system. 

In the case of the deuteron its three e.m. FFs have been 
separated experimentally [35], and it has been possible 
to determine the empirical charge densities. A pictorial 
result for the transverse charge density with transverse 
deuteron polarization sj_ = 1 is shown in Fig. [5] The 
quadrupole field pattern clearly displays a deformation 
along the axis of the spin (cc-axis) together with a small 
overall shift of the charge distribution along the y-axis. 



B. Measuring the shape of hadrons 

After the discussion of these two extreme cases, 
namely, on the one hand of a spin-1 point particle within 
relativistic quantum field theory, and on the other hand 
of a non-relativistic two-body system, we now turn our 
discussion to hadrons, such as mesons and baryons com- 
posed of light quarks. 

Experimentally, accessing information that reveals 
hadron shape, even at the very rudimentary level that 
attempts only to check deviations from spherical sym- 
metry, has proved very difficult for a number of reasons. 
There is only one stable hadron, the proton, and for this 
reason it is the only hadron that can provide a thick 
target for high luminosity precision measurements. The 
relatively long lived neutron either free or inside nuclei 
could provide a possible, but technologically far more dif- 
ficult alternative. Its shape has not been explored so far. 
Both the proton and the neutron are unfortunately spin 
1/2 systems and therefore cannot provide information 
about their intrinsic shape through the measurement of 
a static quadrupole or higher multipole moments. From 
the decuplet spin 3/2 baryons only in the case of the 
A and the it is possible, in principle, to measure 



their quadrupole moments or the transition quadrupole 
moments to some other state. The A + (1232) offers the 
most accessible case; however its exceedingly short life- 
time prevents a viable, yet, experimental way to access 
its quadrupole moment. Nevertheless the magnetic mo- 
ments of the A + [33] and A ++ [30] have been measured, 
albeit with very large errors. New experiments at MAMI 
are expected to yield a more precise measurement for 
the A + dipole moment. The dipole moment of the Sl~ 
is more precisely measured and provides a benchmark 
for lattice QCD calculations 34J, which in turn can pre- 
dict its quadrupole moment. Vector mesons have a static 
quadrupole moment, which, if different from its natural 
value of -1, is a clear indication of a deviation from spher- 
ical symmetry. The p-meson is the lowest lying spin-1 
resonance to test the deviation from spherical symme- 
try. However, experimentally it is again not feasible to 
measure. A beautiful example of what information lat- 
tice QCD can yield on hadron shapes is given in Fig. |6j 
which shows lattice calculations for the density-density 
correlator of the p- meson in the lab frame [35]. In the 
spin projection zero case, the p-meson displays a pro- 
late (cigar-like) deformation in its rest frame. This con- 
clusion is corroborated by a calculation of the p-meson 
quadrupole moment in quenched lattice QCD [36] . 




FIG. 6: Three-dimensional contour plot of the p-meson, of 
spin projection s z =0, density-density correlator (red sur- 
face), showing all positions where the correlator is reduced to 
half its value at the origin. As can be seen, by comparing to a 
sphere (green transparent surface) of radius of approximately 
0.5 fm, the p— meson surface extends outside the sphere at 
the poles whereas at the equator is inside the sphere, showing 
the cigar-like p-meson shape. 

Thus to measure the shape of hadrons, none of the 
"standard" and tested methods used in atomic and nu- 
clear physics can be employed. The only viable path to 
study the nucleon shape remains the one originally pro- 
posed by Isgur and Karl, i.e. to measure the presence of 
resonant quadrupole admixtures in the — > A tran- 
sition, which amounts to determining the off-diagonal 



(transition) quadrupole moment. The theoretical frame- 
work of interpreting these measurements has matured in 
recent years, as will be reviewed below. The precision 
measurements of this transition provide the most reliable 
information we have today for deviation from spherical 
shape for the proton and/or the A + (1232) [3]. 

The experimental technique employed in the determi- 
nation of the deviation from sphericity in the study of 
the de-excitation of the A + (1232) resonance is different 
than those discussed earlier. It involves the detection of 
the radiation pattern of the emitted radiation in the de- 
excitation of the excited state. The concept behind the 
technique derives from classical electromagnetism. The 
observed radiation pattern, its multipole content to be 
precise, reveals information about the shape of the radi- 
ating antenna. The radiation emitted in the de-excitation 
of the the A + (1232) is primarily in the form of pions but 
a small (0.7%) branch of 7 rays is also present. This tech- 
nique of measuring shape rarely has been used in nuclear 
physics, principally due to the experimental complexity 
it presents. An important exemption is the study of the 
15 N excited states using this technique in an (e, e'7) ex- 
periment, which demonstrated both the feasibility and 
accuracy of this method. The 6.33 MeV J p = 3/2~ ex- 
cited state of 15 N, a J p = l/2~ nucleus, presents a case 
where a transition to it from the ground state with the 
same quantum numbers as the 7* TV — > A can be studied. 
The experiment, where the C2 (Coulomb quadrupole), 
E2 (electric quadrupole) and Ml (magnetic dipole) FFs 
were isolated through the tagging of the decay radia- 
tion [37], offers a clear demonstration of the power of 
the technique. The experimental arrangement used is 
shown in Fig. [7] which clearly portrays the concept of this 
experimental technique. The virtual photon causes the 
excitation of the target nucleus and due to angular mo- 
mentum and parity selection rules only magnetic dipole 
and electric quadrupole transitions are allowed; the decay 
radiation pattern allows to identify each admixture. 

In the excitation spectrum of the nucleon the only iso- 
lated state is the A + (1232), which thus allows us to em- 
ploy the same type of measurement as in the 6.33 MeV 
isolated excited state of 15 iV. The 7* A -> A transi- 
tion from J = l/2toJ = 3/2 with no change in parity 
allows us to observe quadrupole E2 and C2 transition 
moments. It is however a mixed transition, which, in 
addition to the quadrupole amplitudes, involves the Ml 
(spin flip) amplitude that is the dominant one. The pres- 
ence of resonant quadrupole strength signifies deviation 
from sphericity of the proton and/or the A + (1232). Us- 
ing the same experimental technique as in the case of N 
it is possible to isolate and measure the weak but impor- 
tant quadrupole amplitudes in the presence of the domi- 
nant Ml transition. Through the extensive study of the 
N — > A transition, pursued during the last thirty years 
using real or virtual photon probes, an extensive body of 
data has emerged that convincingly demonstrates that 
the quadrupole amplitudes are substantial and far bigger 
than can be accommodated by a 'spherical' proton. 




FIG. 7: The detection of the de-excitation radiation pattern 
from a system allows to isolate the contributing multipoles. 
The isolation of the multipole FFs was achieved for the first 
time in the 12 C(e, e'7) and 15 N(e, e'j) reactions where the E2 
and Ml FFs were isolated through the tagging of the decay 
radiation [37] . The figure shows the radiation pattern and the 
isolated FFs for this transition in 15 N. 




60 120 180 



(deg) 

FIG. 8: The precise range of uncertainty that is allowed 
by the Bates and MAMI data for the olt response of the 
p(e, e'p)ir° reaction, at Q 2 = 0.127 GeV 2 and at the posi- 
tion of the A + (1232) resonance, is shown as a function of the 
cm. angle between the proton and virtual photon 9* q . Cur- 
rent phenomenological models (MAID [38] in black, Sato-Lee 
(SL) [39] in blue and DMT gO] in green) predict this sat- 
isfactorily within a 2a confidence level. The corresponding 
calculations for "spherical" nucleon and A (dashed curves us- 
ing same color coding for the various models) cannot describe 
the data; they are excluded to 2 a confidence level. 



Experimentally, the measurement of the small value 
of the electric or Coulomb quadrupole multipole be- 
comes possible through its interference with the dom- 
inant magnetic dipole transition. This is shown in 
Fig. [8] which depicts the angular dependence of the 
longitudinal-transverse interference cross section (<7lt) 
for the p(e,e'p)TT° reaction on top of the A + (1232) reso- 
nance. The cross section a^x is overwhelmingly driven by 
the interference of the dominant transverse Ml amplitude 
with the longitudinal C2 amplitude. The experimentally 



9 



constrained region, by the Bates and MAMI data, is com- 
pared with phenomenological model predictions that at- 
tempt to describe the experimental data. It is evident 
that the model predictions (dotted curves) with reso- 
nant quadrupole amplitudes set to zero, which amounts 
to spherical solutions, are excluded with high confidence. 
The "deformed" model predictions, assuming negligible 
model error, are in agreement at the 2a level, with the 
empirical results. This comparison demonstrates that 
compelling experimental evidence nowadays exists sup- 
porting the conjecture of deformed hadrons. In partic- 
ular, the above data demonstrate with very high confi- 
dence that spherical symmetry for both the nucleon and 
the A + (1232) is experimentally excluded. The experi- 
mental results for hadron deformation in the 7* TV —> A 
transition will be discussed in more detail in Section ITTT1 



C. Calculating the shape of hadrons : lattice QCD 

Having seen clear experimental evidence for a non- 
spherical charge distribution in the N — ¥ A transition, 
we next examine whether this can be calculated and un- 
derstood from QCD, the underlying theory of strong in- 
teractions. 

QCD requires a new methodology in order to evaluate 
quantities related to hadron structure, the reason being 
that hadrons are bound state systems having a mass that 
is mostly generated by the interaction rather than by the 
sum of the mass of their constituents. Perturbative QCD 
has been very successful in describing high energy pro- 
cesses. On the other end of very low energy, chiral pertur- 
bation theory has provided the appropriate effective field 
theory framework for precise calculations of observables 
in terms of a small expansion parameter, such as an exter- 
nal momentum or pion mass. This framework provides 
a systematic expansion involving an increasingly large 
number of low-energy constants (LECs). The latter are 
free parameters, which are beyond the predictive power 
of the effective field theory. Some of these have been de- 
termined from phenomenological information, however, 
the vast majority remains unknown limiting the predic- 
tive power of chiral effective field theory. To calculate 
LECs from the underlying theory of QCD as well as to 
make predictions beyond a regime where a perturbative 
or small scale expansion is applicable requires an inher- 
ently non-perturbative technique. Such an approach that 
enables us to solve the theory in the non-perturbative 
domain starting from the underlying QCD Lagrangian 
is lattice QCD, a discretized version of QCD formulated 
in terms of Feynman's path integrals on a space-time 
lattice preserving gauge symmetry |41j . Like the con- 
tinuum theory, the only parameters are the bare quark 
masses and the coupling constant. One recovers contin- 
uum physics by extrapolating results obtained at finite 
lattice spacing a to a = 0. 

A crucial step, that enables one to numerically eval- 
uate the path integrals needed, is rotation to imaginary 



time, t — > —it, resulting in the replacement of the time 
evolution operator exp(—iHt/h) by exp(—'Ht/h). Within 
the Feynman path integral formulation, observables are 
calculated by a weighted sum over all possible trajecto- 
ries. In imaginary time it becomes possible to generate a 
representative ensemble of trajectories by using stochas- 
tic methods analogous to those applied in the evaluation 
of observables in statistical mechanics. 

Calculations in lattice QCD started in the early 80's, 
and during the first two decades were performed mostly 
in the quenched approximation, which neglects pair cre- 
ation. This enormously simplifies the generation of the 
gauge fields via Monte Carlo methods since one is left 
with a local gauge action. During the past ten years, the- 
oretical progress in combination with terascale computers 
have made simulation of the full theory with light pions 
and large enough volumes feasible using several different 
discretization schemes. The simplest lattice QCD action 
is due to Wilson [41]. Nowadays, one uses improved dis- 
cretized versions of the Dirac operator with reduced finite 
lattice spacing artifacts and better chiral properties, all 
of which are expected to yield the same results in the 
continuum limit |42) . Using these improved fermion dis- 
cretization schemes, simulations with pion masses within 
100 MeV of the physical pion mass are currently available 
with simulations using improved Wilson fermions even 
reaching the physical pion mass [43 . A benchmark cal- 
culation for lattice QCD is the evaluation of the low lying 
hadron spectrum, where a systematic study of the hadron 
masses using different discretization schemes has been 
performed and the continuum and infinite volume limits 
have been examined. The agreement with experiment ob- 
served from such systematic lattice studies [121 [H] , pro- 
vides a validation of the lattice QCD approach, paving 
the way to use lattice QCD to provide predictions for 
quantities, which are very difficult to access experimen- 
tally, as for example the e.m. FFs of an excited hadronic 
state, such as the A(1232) resonance. It furthermore al- 
lows to study how the physics is affected when varying 
fundamental parameters such as quark masses outside 
their values realized in nature. 

Information on hadron shapes can be extracted from 
FFs and generalized parton distributions. The evaluation 
of these quantities is more involved than the computation 
of hadron masses. FFs are connected to hadron matrix 
elements of the type (h! (p')\0\h(p)) and one needs, in 
general, to compute the diagrams shown in Fig. [9j where 
the solid lines denote fully dressed quark propagators. 
The diagram where the operator couples to a sea quark, 
shown in upper panel of Fig. [9j is particularly difficult to 
calculate since it involves a disconnected quark loop. For 
the evaluation of transition FFs where the final hadron 
state h! has different quantum numbers from the initial 
h, the disconnected diagram vanishes. For diagonal ma- 
trix elements, assuming isospin symmetry, the discon- 
nected contribution vanishes for isovector operators and 
therefore isovector FFs can be calculated from the con- 
nected diagram alone. Although recently efforts to cal- 
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FIG. 9: Upper panel: Disconnected; Lower panel: connected 
diagrams. With h and h' we denote hadronic states and with 
O the operator of interest. 



culate such disconnected contributions have intensified, 
up to now lattice computations of FFs generally neglect 
disconnected contributions. The standard procedure to 
evaluate the connected three-point function shown in the 
lower panel of Fig.[9]is to compute the so called sequential 
propagator, a convolution of fixed source quark propaga- 
tors, which are technically straight forward to calculate. 
In most recent studies of the e.m. FFs the fixed sink 
method is used. Its name comes from the fact that we 
consider a given final state h' is at a fixed time £2 from 
the initial state h created at time zero. Within this ap- 
proach, any operator can be inserted at any intermediate 
time slice ti, as seen in Fig.[9j carrying any possible value 
of the lattice momentum. For a recent review, which in- 
cludes comparison of nucleon electromagnetic FFs within 
a number of different discretization schemes see Ref . |45j . 

To probe hadron deformation, the e.m. current is used 
for the operator O in Fig.[9| Lower moments of transverse 
spin densities of quarks in the nucleon |46| or pion [47] 
as well as the transverse momentum dependent parton 
distribution functions |45j can also be evaluated using 
these techniques. 

As an example of the predictive power of lattice QCD, 
we show in Fig. 10 the transverse charge density of Eq. ^ 
for a A + (1232jthat has a transverse spin projection 
,s_l = +3/2. This charge density is obtained from the 
A e.m. FFs, calculated within lattice QCD in Ref. [48] . 
as described in more detail in Section HVl As can be seen 



from Fig. 10 the quark charge density in a A + in a state 
of transverse spin projection sj_ = +3/2 is elongated 
along the axis of the spin (prolate deformation) when 
observed from a light-front. 

Although lattice QCD provides an ab initial calcula- 
tion of fundamental quantities such as FFs or moments of 
generalized parton distributions, a careful analysis of sta- 
tistical and systematic errors must be performed before 
one can reliably compare to experiment. The systematic 
errors arise because lattice calculations necessarily are 



FIG. 10: Lattice QCD results for the quark transverse charge 

density p£j 3 in a A + (1232) which is polarized along the posi- 
1 2 

tive x-axis. The light (dark) regions correspond to the largest 
(smallest) values of the density. In order to see the defor- 
mation more clearly, a circle of radius 0.5 fm is drawn for 
comparison. The density is obtained from quenched lattice 
QCD results at m n = 410 MeV for the A e.m. FFs [48]. 



performed for a finite lattice size and spacing a due to 
the discretization of space-time, which breaks continuous 
rotational invariance to a discrete one. These systematic 
errors need to be investigated by repeating the calcula- 
tion for various volumes and lattice spacings. Except for 
hadron masses, where both the infinite volume and zero 
lattice spacing limits are taken, for other quantities like 
FFs such an analysis has just began. Another source 
of systematic error is the fact that FFs calculations still 
utilize dynamical quarks of larger mass than the physical 
one. Whereas finite a, lattice size L and magnitude of the 
quark masses are amenable to systematic improvements, 
rotation to Euclidean space selects a set of observables, 
determined from the properties of the discrete low-lying 
states, which can be studied within this framework. Nu- 
cleon FFs and moments of parton distributions are exam- 
ples of such observables. Excited states are more difficult 
to compute since they are exponentially suppressed as 
compared to the ground state due to the Euclidean time 
evolution. Techniques have been developed to extract 
the low lying excited states, however most calculations 
are still done in quenched QCD and without an analy- 
sis of systematic errors, although some recent results on 
the excited states of the nucleon using two dynamical 
quarks have been presented. The study of resonances in 
lattice QCD is a recent activity. One of the reasons is 
that up to very recently the quark masses that could be 
simulated were too large to allow decays. Although ex- 
traction of the spectral function from lattice correlators 
is not feasible since the low energy continuum scattering 
states dominate, there are theoretical techniques to study 
the width of resonances [49 that make use of the depen- 
dence of the energy on the finite lattice length. These 
techniques, combined with the background field method, 
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can yield the magnetic and electric quadrupole moments 
of resonant states [50]. However, to go beyond the calcu- 
lation of the decay width and the lower moments to the 
calculation of FFs for resonances such as the A is still an 
open theoretical problem. 

In Scction |TV} we will present results showing the state- 
of-the-art of the lattice calculations for the e.m. FFs of 
the A(1232) resonance, as well as for the e.m. FFs de- 
scribing the 7* TV — > A transition, and discuss the result- 
ing theoretical predictions for hadron deformation. 

III. EXPERIMENTAL EVIDENCE 

The experimental landscape concerning the investiga- 
tion of the shape of hadrons has been dominated by 
the quest for resonant quadrupole amplitudes in the 
7* A — > A transition in the proton. Recently other re- 
actions have been suggested, e.g. the study of the 
7* A — ¥ A transition in neutrons or in nuclei, and they 
may become technically feasible in the near future. In 
addition, it is understood that the detailed and precise 
understanding of form factors can bring new complemen- 
tary information on the issue of the shape of hadrons. In 
addition to the formidable technical difficulties of access- 
ing new reaction channels to the required precision, the 
theoretical framework to extract the important physical 
conclusions needs to be further developed. 

The experimental investigation of the 7* N —> A tran- 
sition can be classified according to the reaction channel 
probed. The A + (1232) can be excited by real or virtual 
photons, 7* , and decays through pion or photon emission: 

7 *p A+(1232) ^ptt° (66%), 
7 *p ->■ A+(1232) mr+ (33%), 
7*p ->■ A+(1232)^p 7 (0.56%). 

The pion decay channels have been extensively explored 
while the third, involving the gamma decay branch, has 
been studied with real compton scattering (RCS). Virtual 
compton scattering measurements (VCS) are beginning 
to emerge with the aim of mapping the polarizabilities 
at high missing mass |51j and/or investigating the issue 
of deformation [S"2"] . 

The first generation 7* A — > A and in general nucleon 
resonance experiments were conducted in the late six- 
ties and early seventies, before the issue of deforma- 
tion was even raised, at DESY, NINA and CEA with 
low quality beams and experimental equipment not de- 
signed to address such refined questions. The data 
that emerged were characterized by limited accuracy, 
but they did provide valuable guidance on the design 
of future experiments [53] . The second generation ex- 
periments were obtained by a newer generation of ac- 
celerators at Brookhaven, Bates, MAMI, and CEBAF 
with optimized equipment and in general with polarized 
beams. Third generation experiments are now beginning 
to emerge; they have been conducted primarily with po- 
larized and tagged real photons, impinging on polarized 



targets. Electroproduction experiments with polarized 
targets are particularly difficult with only one measure- 
ment reported in the literature using the internal target 
facility at NIKHEF [M] and having low statistical accu- 
racy. The JLab Hall A experiment |55j . which presented 
high quality extensive recoil polarization measurements 
using polarized beams, is a truly third generation ex- 
periment, which both demonstrated the feasibility of the 
technique, the precision that can be achieved and the rich 
physics output that can emerge. 

In general, in the real photon sector, the "second gener- 
ation" experiments arc completed and analyzed and the 
era of "third generation" experiments is about to begin in 
earnest, in view of the important instrumentation initia- 
tives [5(1] at Mainz and at Bonn. JLab and MAMI C have 
optimal beams and detection systems for the pursuit of 
this program, which is far from being exhausted. 

The "deformation" signal in the real photon sec- 
tor comes from the study of the transverse electric 
3 li 

quadrupole (E 1+ also denoted by E2) multipole. If vir- 
tual photons are used, the longitudinal quadrupole (L\^ 
or C2) is also accessed. The superscript indicates the to- 
tal isospin 3/2, whereas the subscript denotes the I = 1 
angular momentum in the ttN system, and the "+" refers 
to the total angular momentum J = I + 1/2 = 3/2. In a 
quark-model picture, the 7-/V — > A transition is described 
by a spin flip of a quark in an S-wave state in the nu- 
cleon, resulting in a magnetic dipole (M^{ 2 or Ml) tran- 
sition. Any D-w&ve admixture in the nucleon or the A 
wave functions, also allows non-zero values for the elec- 
tric quadrupole (C2 and E2) transition. This is depicted 
graphically in Fig. [TTJ 

P(938) A(1232) 




M 1+ ,E 1+ ,S 1+ n' 



Spherical => M1 Deformed => M1 , E2 , C2 

FIG. 11: Quark model picture of Ml, E2 and C2 amplitudes 
in the N — > A transition induced by the interaction of a 
photon (real or virtual) with a single quark in the nucleon. 
Presence of quadrupole amplitudes in the transition requires 
N and/or A wave functions to have a D-wave component 
(indicated by a non-spherical shape). 

It has become standard practice in the field to mea- 
sure the resonant quadrupole strengths relatively to 
the resonant dipole by introducing the ratios EMR= 

ImE^/JxaM^ and CMR= ImL^ 2 /ImAffj 2 . EMR 
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— energy dep (global) fit, Hanslein etal 
* energy indep, (local) fit. Hanslein el ai 
t exp analysis. Krahn etal, Mainz 1997 
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FIG. 12: Results from MAMI (left) and LEGS (right) have yielded precise measurements of the resonant quadrupole amplitude 
at the photon point. A most sensitive probe is the polarization asymmetry E which has been measured precisely at both MAMI 
and LEGS. The derived multipoles from the MAMI cross sections, yield an accurate measurement of EMR. 



and CMR have thus become the signal of deformation. 

A. Real Photon Measurements 

In photoproduction, the presence of a resonant 
quadrupole amplitude, is particularly hard to isolate be- 
cause the transverse channel is overwhelmed by the mag- 
netic dipole (Ml) amplitudes and contaminated with 
other "background" (non resonant) processes of simi- 
lar magnitude. In this sense, the appears in next 
to leading order (NLO) in photoproduction. Precision 
measurements with polarized tagged photons performed 
at Mainz (MAMI) and Brookhaven (LEGS) in the late 
nineties represent a tour de force of experimental fi- 
nesse. The small quadrupole amplitude has been de- 
tected in the measurement of the polarization asymme- 
try E = (crij — tTj_ ) / (cri + ax) shown in the right panel of 
Fig. |12| The asymmetry E is measured with reduced sys- 
tematic error by flipping the polarization of the tagged 
photon beam parallel (||) and perpendicular (_L) to the 
scattering plane. Analysis of the MAMI (j,tt + ) and 
(7, 7T°) data, yields the impressive results shown in the 
left panel of Fig. [l2j It is obvious from this figure that 

the derived results heavily depend on the W dependence 

3/2 

of the cross section. The E 1+ multipoles have a striking 
non-resonant shape, a manifestation of the complicated 
processes that contribute to this channel. The measure- 
ments from MAMI [57 and LEGS [58 have converged 
as far as the determination of the asymmetries are con- 
cerned. The resulting EMR values are : 

LEGS : EMR = -(3.07 ± 0.26 sta t. +sys t. ± 0.24 mod .)%, 
MAMI : EMR=-(2.5±0.W±0.2 syst .)%. 



A number of theoretical calculations are in good agree- 
ment with the experimentally derived EMR value. Both 
the (7, 7T°) and the (7,71+) channels have been studied 
extensively. The (7, 7) channel (RCS) has also been stud- 
ied [59] . where the resonance pion-photoproduction am- 
plitudes were evaluated leading to the multipole EMR 
(340 MeV) =(-1.6±0A {stat+syst) ±0.2 (model) )%, in rea- 
sonable agreement with the photopion measurements. 

The situation concerning the — > A transition in the 
real photon sector has remained stable, without experi- 
mental results reported to change this picture in the last 
five years. A subsequent analysis [60] and new data |61j 
give EMR = (-2.74±0.03- tQi )±0.3 (syst) )%, confirming 
the EMR values of Ref. [5?, 58 . In the closely related 
areas of threshold pion production [6 2) and in the mea- 
surement of the magnetic dipole [33, 63J of the A + (1232), 
the precise results that emerged provide both a test as 
well as valuable guidance to theory and phenomenology 
that is common to both. The installation of the Crystal 
Ball at MAMI and of the Crystal Barrel at ELSA, have 
brought new very powerful tools, which are expected to 
yield even more precise data and results. 

B. Electroproduction measurements 

In electron scattering experiments, in addition to the 
transverse responses, the longitudinal responses are also 

accessible, which are sensitive to leading order to the 

3/2 

longitudinal quadrupole multipole, L 1+ or C2. Further- 
more, the Q 2 evolution of the various responses offers the 
ability to distinguish between different processes. This is 
of particular value for the understanding of the distinc- 
tive roles played by the mesonic cloud as compared with 
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the quark core. These advantages are technically chal- 
lenging and time consuming to realize due, primarily, to 
the numerous measurements needed to cover the widest 
possible range of momentum transfers. 

Consistent results have been reported from several 
groups [SSJ IMH7T] at Bates, ELSA, MAMI and JLab 
mapping the momentum transfer range from Q 2 — 0.06 
to 6.0 GeV 2 with high precision in a limited number of 
observables sensitive to the issue of deformation. How- 
ever, there are still discrepancies on the extracted EMR 
and CMR values, which are not directly measurable due 
to the methodology used in extracting multipoles, an is- 
sue discussed in the next section. 

Starting from the experimental observables, two meth- 
ods have been used for extracting multipole amplitudes: 
a) The Truncated Multipole Expansion (TME) approx- 
imation in which most or all of the non-resonant mul- 
tipoles are neglected (e.g. see [Ml [SB]) assuming that, 
at resonance, only the resonant terms contribute signifi- 
cantly and are fitted to the data, and b) the Model De- 
pendent Extraction (MDE) method where a phenomeno- 
logical reaction framework with adjustable quadrupole 
amplitudes is used, e.g. see [66l [67l [73]. It is assumed 
that the reaction is controlled at the level of precision 
required for the disentanglement of the background from 
the resonance amplitudes. Clearly the MDE method is 
superior, given the sophistication that phcnomenological 
models have achieved in describing the data. 

In the recent electroproduction experiments, which al- 
most invariably are carried out with polarized beams, 
the transverse-longitudinal response functions cttl an d 
cttl' are measured. Their simultaneous measurement al- 
lows the extraction of the real and imaginary parts of 
the same combination of multipole amplitudes. Knowl- 
edge of both responses is particularly valuable because 
<ttl is most sensitive to the presence of a resonant longi- 
tudinal quadrupole amplitude, while cttl' is particularly 
sensitive to the background contributions, thus providing 
information on the two aspects of the problem that need 
to be controlled independently [75]. The importance of 
background is clearly seen in the W behavior of the re- 
sponses [67] and the non-vanishing recoil polarization 
P n [551 HFJ j which bears close resemblance to cttl' • The 
transverse-transverse response, cttt, which is sensitive to 
the electric quadrupole amplitude, was only recently iso- 
lated for the first time at non-zero Q 2 , with experiments 
pursued at Bates, JLab and MAMI [751150] . 

Fig. [13] offers a compilation of CMR and EMR as a 
function of Q 2 . Both EMR and CMR are small and neg- 
ative in the region where they have been measured. From 
the accuracy of the present data one immediately recog- 
nizes that the quark model predictions, which historically 
provided the motivation for these measurements, do not 
agree with the data. In particular the dominant Ml ma- 
trix element is found to be ~ 30% stronger and the E2 
and C2 amplitudes at least an order of magnitude larger 
and often of a different sign than the predictions of quark 
models. This failure is to be expected since the quark 
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FIG. 13: The experimentally derived values for M^ 2 (Ml) 
(top panel), CMR (middle panel) and EMR (lower panel) 
compared to phenomenological model results for Ml, and to 
lattice for EMR and CMR. The derived multipole ratios are 
shown without the model error that is of the order or larger 
than the depicted experimental error. 



model does not respect chiral symmetry whose dynamic 
breaking leads to a strong, non-spherical, pion cloud sur- 
rounding hadrons [81) . As will be discussed in detail in 
Section [TV] it was realized that the pion cloud was a 
necessary ingredient to be added to quark models. This 
is demonstrated by the calculations of the SL [39] and 
DMT [?0] models, which describe the data adequately, 
and which show that most of the strength of the re- 
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sponses (and the EMR and CMR values) at very low 
Q 2 values, below ~ 0.25 GeV 2 /c 2 , arises on account of 
the mesonic degrees of freedom. The recent results from 
MAMI along with the earlier ones from Bates [75] and 
the recent low Q 2 measurements from CLAS [5U], give 
strong support to this interpretation. At asymptotic val- 
ues of Q 2 helicity conservation 82J requires that EMR 
— > 1 and that CMR — > constant. Clearly this regime has 
not been reached. The upgrade of CEBAF to 12 GeV 
will allow to extend the measurements to higher Q 2 , al- 
though this will pose significant challenges in isolating 
the relevant partial cross sections and even bigger ones 
in extracting the relevant amplitudes. 

Finally, the H (e, e'p) 7 channel (VCS), which only re- 
cently has been accessed, allows the extraction of the 
quadrupole amplitudes through a purely electromagnetic 
reaction channel providing an important cross check to 
the derived results from the pionic channel. The disper- 
sion theory of Ref. [53] , allows one to address the physics 
of deformation and of nucleon polarizabilities in the re- 
gion above pion threshold simultaneously. Recent results 
from MAMI report the extraction of polarizabilities [3T] 
and the first observation of VCS data sensitive to the 
resonant quadrupole amplitudes [52] ■ The results are 
in excellent agreement with those derived from the pion 
channel. 



C. Sensitivity, Precision and Estimation of 
Uncertainties 

The 7* A A data up to Q' 2 = 6.0 GeV 2 are, in gen- 
eral, characterized by small systematic errors and high 
statistical precision. The interpretation of the data in 
terms of the deformation has been demonstrated, and as 
a result, the research thrust shifted from the investiga- 
tion of whether the conjecture for deformation is valid 
to the exploration of the mechanisms that cause it. In- 
vestigating the physical origin of deformation requires 
the measurement of new responses and the comparison 
of the theoretical results with the experimentally derived 
quantities, at a level of precision far superior to the one 
feasible today. This detailed comparison necessitates a 
reliable determination of the uncertainties of both the 
experimental results and the theoretical calculations. 

The need for a critical and precise comparison of data 
and theory when extracting multipoles in nucleon reso- 
nance studies is reminiscent of the "crisis" in the analy- 
sis of electron scattering data in the early 1970s, where 
the very precise data could not be meaningfully com- 
pared with the theoretical calculations in order to de- 
rive nuclear charge densities. This was primarily due to 
the lack of an appropriate methodology that could en- 
able to quantify the uncertainties in the extracted densi- 
ties, which, like multipole amplitudes, are not experimen- 
tal observables. The resolution of the "crisis" through 
the introduction of a "Model Independent" extraction of 
charge densities led to a revolution in the field and to the 



outstanding achievements in electron scattering. 

The leading method of extraction of multipole am- 
plitudes, the Model Dependent Extraction (MDE), pro- 
duces extracted values that are biased by the model and 
characterized by a model error, which is hard to estimate, 
especially if a single model is employed [B"6"l I6"5l 1701 [721 174] . 
An Ansatz for estimating the model uncertainties in the 
extracted multipoles has been proposed [84] [85] and used 
in a few cases [7TJ [73] . In this method the same data 
are analyzed employing different models, which describe 
the data adequately, and attributing the resulting spread 
in the extracted quantities to model uncertainty. Even 
though the phenomenological models available are of con- 
siderable sophistication, the small non-resonant ampli- 
tudes collectively could induce large correlations and er- 
ror in the extraction of the resonant amplitudes, result- 
ing in the unsatisfactory situation that the uncertainty 
is only approximately known. 

A novel model independent method, the Athens Model 
Independent Analysis Scheme (AMIAS), for extracting 
multipole information from experimental nucleon reso- 
nance [55] and for analyzing lattice QCD simulation data 
has been presented [57]. The method quantifies the un- 
certainty of the extracted multipoles and yields new in- 
formation on background amplitudes, which MDE is in- 
capable of accessing [55j . Results from AMIAS are shown 
in Fig. [8] (yellow bands) for the CMR sensitive <jlt par- 
tial cross section with the precisely defined one-er uncer- 
tainty. The experimentally allowed o~lt partial cross 
section, as constrained by the Bates and MAMI data at 
Q 2 = 0.127 GeV 2 , are shown as a function of 6* . They 
are compared with theoretical model predictions that ac- 
count for them. It is evident that the model predictions 
(dotted curves) with resonant quadrupole amplitudes set 
to zero, which amounts to spherical solutions, are ex- 
cluded with high confidence. On the contrary model 
predictions from models that allow mesonic degrees of 
freedom allowing for deformation are in reasonable agree- 
ment with the experimental results at the 2a level. Dif- 
ferences among the curves predicted by the various phe- 
nomenological models are visible, but no inference can be 
drawn as their model error is not known. Nevertheless, 
the comparison demonstrates with extremely high con- 
fidence, with experimental errors precisely defined, that 
the assumption of sphericity for both the nucleon and the 
A + (1232) is incompatible with the data. 



IV. THE SHAPE OF NUCLEON AND A 
RESONANCE : THEORETICAL 
UNDERSTANDING 

Having seen first experimental evidence for a non- 
spherical shape of the nucleon and A resonance, we next 
discuss its theoretical understanding. For the A reso- 
nance, its e.m. FFs are not accessible experimentally 
They are, therefore, an ideal example of observables 
where lattice QCD can make predictions. The state-of- 
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the-art of these calculations as well as their implication 
on the shape of the A resonance are discussed. Sub- 
sequently our current theoretical understanding of the 
7* A — > A transition is summarised including an inter- 



pretation of the data presented in Section III 



A. A charge densities : lattice QCD 

As the nucleon is a spin-1/2 particle, its transverse 
charge densities do not exhibit a quadrupole pattern, nor 
do they encode any information on its shape. For spin- 
3/2 baryons, such information can however be obtained 
from the charge densities. 

The matrix element of the e.m. current operator 
J M between spin-3/2 states, such as the A(1232) res- 
onance, can be decomposed into four multipole transi- 
tions: Coulomb monopole (EO), magnetic dipole (Ml), 
Coulomb quadrupole (E2) and magnetic octupole (M3), 
described by the corresponding FFs Geo, Gmi, Gei and 
G M 3 [13 W\- Thcir values at Q 2 = define e.g. the 
magnetic dipole moment : /ia = Gj\/i(0)e/(2MA), or 
the electric quadrupole moment : Qa = G£ 2 (0)e/M A . 

The empirical knowledge of the A electromagnetic mo- 
ments is scarce, even though there were several attempts 
to measure its magnetic dipole moment. The current 
PDG value of the A + magnetic dipole moment is [50] : 

H A + = 2.7+J-°(stat.) ± 1.5(syst.) ± 3(theor.) /Mr . (13) 

This result was obtained from radiative photoproduc- 
tion (7 A — > irN'y') of neutral pions in the A (1232) region 
by the TAPS Collaboration at MAMI [53], using a phe- 
nomenologica l m odel of the jp — > ir°pj' reaction. For 
the A+, Eq. fl3]) implies : 



Gmi(O) = 3.5I^(stat.)±2.0(syst.)±3.9(theor.) . (14) 

The size of the error-bar is rather large due to both ex- 
perimental and theoretical uncertainties. 

For the A electric quadrupole moment or magnetic oc- 
tupole moments, no direct measurements exist, nor do 
we have any empirical information on the Q 2 behavior 
of the A e.m. FFs. We thus rely on recent lattice QCD 
calculations [48] that can predict these FFs. 

Calculation of the e.m. FFs within lattice QCD re- 
quires the evaluation of a three-point function, as de- 
picted in Fig. [9] We only consider here the connected 
diagram. Its evaluation involves two spatial sums : one 
over the spatial coordinates of the operator and one over 
the spatial coordinates of the final state. In the so-called 
fixed sink method, the sum over X2 is done automatically 
by generating a sequential (backward) propagator from 
the sink to the operator. Inserting the operator, which 
can be done at all values of x\, and summing over x\ with 
the appropriate Fourier phase and propagator starting at 
t = and ending at t = t\ yields the connected three- 
point function, for all momentum transfers q. To extract 
the matrix element (h' (p')\0\h(p)) one studies the large 



t\ Euclidean time behavior of an appropriately defined 
ratio of the three-point function and two-point functions, 
yielding a time independent quantity (plateau). Such a 
behavior signals identification of the lowest hadron states 
h and h! from the tower of QCD states with the same 
quantum numbers as h and hi . Fitting to this plateau 
value we extract the matrix element (h' (p')\0\h(p)) and 
from this, depending on the choice of O, the FFs or mo- 
ments of parton distributions. 
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FIG. 14: Lattice results for the A + (1232) form factor Gmi 
(upper panel) and Gei (lower panel) at the smallest pion 
mass in three simulations [48]. The lines show the fits to an 
exponential form of the quenched lattice results and to the 
results obtained using the hybrid action. The error band is 
calculated using a jackknife analysis on the fitted parameters. 

To probe hadron deformation we will use in the fol- 
lowing the e.m. current as the operator O. Since the 
connected diagram for the A e.m. FFs is calculated by 
performing sequential inversions through the sink, the 
initial and final A states need to be fixed. The A is de- 
scribed by a Rarita-Schwinger spinor and therefore there 
is some freedom in the vector indices that can be cho- 
sen. Ref. [48 concentrates on a few carefully chosen 
combinations that best determine the three FFs pay- 
ing particular attention in constructing a combination 
that isolates the electric quadrupole FF. In order to ef- 
ficiently check the lattice set-up, a quenched calculation 
is carried out using Wilson fermions and the standard 
Wilson plaquette gauge action [35] for which statistical 
fluctuations are small. Quenched results are then com- 
pared to a calculation using two dynamical degenerate 
flavors of Wilson fermions (Np = 2) and the standard 
Wilson plaquette gauge-action as well as using a hybrid 
action 48 . The latter case uses two degenerate flavors 
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of light staggered sea quarks and a strange staggered sea 
quark (Np = 2 + 1) simulated using the Asqtad MILC 
action [5D] . The strange quark mass is fixed to its physi- 
cal value. These gauge configurations are among the best 
simulations of the QCD vacuum available. The valence 
quarks are domain wall fermions (DWF) that preserve a 
form of chiral symmetry on the lattice. A comparison 
between results obtained with these two different lattice 
formulations for the quarks (i.e. dynamical Wilson and 
staggered sea with DWF) provides a non-trivial check of 
lattice artifacts. In both dynamical simulations the A is 
a stable particle. 



We show the results for the A FFs Gmi and Gei in 
For the pion masses considered, there is agree- 



Fig. 14 



ment among results using the different actions, with sta- 
tistical errors being smallest in the quenched theory, as 
expected. For the A magnetic dipole moment, first dy- 
namical results, using a background field method, with 
Np = 2 + 1 quark flavors were presented in Ref . [50] . The 
magnetic moment can also be extracted by fitting the (in- 
dependence of the magnetic dipole form factor Gmi to 
determine its value at Q 2 = 0. The values obtained in 
these two approaches are in agreement |48j . 



Having a determination of the A e.m. FFs in lattice 
QCD one can calculate its transvers e lig ht-front charge 
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Choosing the 



density p^ s± [48 , as shown in Fig. 
transverse spin vector S± = e x , the electric quadrupole 
moment in a state of sj_ = +3/2 for such charge distri- 
bution is then obtained from Eq. ^ as : 



Qa 



2 [2 (Gmi (0) 



3e A ) + (G £2 (0) + 3e A )] 



Ml 



(15) 

Note that for a spin-3/2 particle without internal struc- 
ture, for which Gmi(0) = 3e A and Ge 2 (0) — ~ 3e A , the 
quadrupole moment of its transverse light-front charge 
density vanishes. This is in contrast with the non- 
relativistic case, where a non-zero value of Gei is usually 
interpreted as a non-zero quadrupole moment in the lab 
frame. It is thus interesting to observe from Eq. ( 15 ) 
that, as for the case of a spin-1 particle discussed in Sec- 
tion II, Q s± is only sensitive to the anomalous parts of 
the spin-3/2 magnetic dipole and electric quadrupole mo- 
ments, and vanishes for a particle without internal struc- 
ture. Extrapolating the A e.m. FFs to Q 2 = 0, and using 
the extracted values in Eq. ( [15] ) yields a quadrupole mo- 
ment <9 + 3, of (0.73 + 0.16) (ejMD for the quenched and 
(0.51 ± 0.22) (e/M A ) for the hybrid cases. Both calcula- 
tions therefore show a (small) prolate deformation of the 
two-dimensional light-front charge density along the axis 
of the A spin (for the case of spin projection +3/2). 



B. The electromagnetic N — > A transition in QCD 

1. Electromagnetic moments and densities 

Direct experimental evidence for a deformation of N 
and A states can be obtained from the 7* TV A transi- 
tion, which is usually characterized in terms of three 
Jones-Scadron FFs [HI] : G* M1 , G* E2 and G* C2 , denoting 
the magnetic dipole, electric quadrupole and Coulomb 
quadrupole transitions respectively. For a review and 
more details see Ref. [53] ■ In the following, we will often 
discuss the ratios EMR and CMR, which are expressed 
in terms of the Jones-Scadron FFs as: 



EMR = - 



G M^ i C MR=-%-^@ 2 -, (16) 
mi 4M| G* M1 



G 



with Q± = y/ (M A ± M N ) 2 + Q 2 . 

From the experimental information on the 7*jVA tran- 
sition, discussed in Section [ill] one can extract the tran- 
sition magnetic dipole and electric quadrupole moments 
from the values of the FFs at O 2 = [52]: 



Hn^a = y/M&/M N G* M1 (0) \p, N ], 

2M A 



Q 



]V->A 



-6y/M A /M N - 



M N {Ml-M 2 

Using the experimental information, this yields 

Mp-s-A+ = [3.46 ± 0.03] /z at, 
Q P ^A+ = -(0.0846 + 0.0033) fm 2 . 



(17) 

72T^ 2 (0)(18) 



(19) 
(20) 



One often uses an equivalent parametrization for the 
7iVA transition at the real photon point (Q 2 = 0) 
through two helicity amplitudes A1/2 and A 3 / 2 , where 
the subscript denotes the total 7 + N helicity in the A 
rest frame. Furthermore, one can generalize the consid- 
erations for the nucleon and A FFs to extract from the 
empirical information on the Q 2 dependence of the Ml, 
E2, and C2 transition FFs, the quark transition charge 
densities in the transverse plane, which induce the e.m. 
N — >• A excitation [53] . The transition density in a trans- 
versely polarized N and A shows both monopole, dipole, 
and quadrupole patterns. The latter, shown in Fig. [15] 
maps the spatial dependence in the deformation of the 
transition charge distribution. 



2. Model descriptions of the 7* TV A transition 

As discussed above, the e.m. TV — > A transition is pre- 
dominantly of the magnetic dipole (Ml) type. A first 
understanding of the 7*./VA transition can be obtained 
based on symmetries of QCD and its large number- 
of-color (N c ) limit. In this limit, the baryon sector 
composed of up, down, and strange quark flavors of 
QCD displays an SU{6) spin- flavor symmetry. This 
spin-flavor global symmetry of QCD is at the basis of 
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FIG. 15: Quadrupolc contribution to the transverse charge 
density for the p — > A + transition, when iV and A are po- 
larized along the x-axis with spin projection +1/2. For the 
N -> A e.m. FFs, the phenomenological MAID2007 [38] 
parametrization is used. Figure from Ref. [23] . 

many quark models, in which baryons are described as 
(non-relativistic) quantum-mechanical three-quark sys- 
tems moving in a confining potential. In such quark- 
model picture, the N — > A transition is described by an 
Ml spin flip of a quark in the S'-wave state, illustrated 
The SU(6) symmetry allows to relate the 



in Fig. 11 



magnetic dipole moments of the proton and the p — > A + 
transition as : /i p ^A+ = 2y/2/3fip — 2.63 /in, which 
is about 25 % lower than the experimental number of 
Eq. (19). Any _D-wave admixture in the nucleoli or the 



A wave functions allows non-zero values for the E2 and 



C2 quadrupole transitions, as illustrated in Fig. 11 

The prototype quark model is the Isgur-Karl 
model |93j . where the constituent quarks move in a 
harmonic oscillator type long-range confining potential, 
which is supplemented by an interquark force correspond- 
ing with one-gluon exchange. The one-gluon exchange 
leads to a color hyperfine interaction, which was found 
to predict well the mass splittings between octet and de- 
cuplet baryons 94 . This hyperfine interaction contains 
a tensor force which produces a D-state admixture in the 
N and A ground states, around 1 % [TSJ |SS] . As a result 
of such D-wave components, the N and A charge densi- 
ties become non-spherical, yielding small negative EMR 
values, in the range -0.8% < EMR < -0.3% within 
non-relativistic quark models [T51 [35] . The small value 
for EMR already indicates that any effect of deformation 
in the nucleon and/or A ground state is rather small and 
very sensitive to details of the wave function, as well as 
truncation in the quark model basis [HT] [§E] ■ The error 
induced due to the truncation in the quark model ba- 
sis has been further investigated in the relativized quark 
model [99| , typically resulting in an even smaller negative 
value, namely EMR ~ -0.2%. 

Even though the constituent quark model, despite its 
simplicity, is relatively successful in predicting the struc- 
ture and spectrum of low-lying baryons, it under-predicts 
/ijv^A by more than 25 % and leads to values for EMR, 
which are typically smaller than experiment. More gen- 



erally, constituent quark models do not satisfy the sym- 
metry properties of the QCD Lagrangian. In the limit 
of massless up and down (current) quarks, the QCD La- 
grangian is invariant under SU(2)l x SU(2)r rotations 
of left (L) and right (R) handed quarks in flavor space. 
This chiral symmetry is spontaneously broken in nature 
leading to the appearance of massless Goldstone modes, 
pions, which acquire a mass due to the explicit breaking 
of chiral symmetry. Since pions are the lightest hadrons, 
they dominate the long-distance behavior of hadron wave 
functions. As the A(1232) resonance nearly entirely de- 
cays into 7tjV, the pions are of particular relevance to the 
~/*NA transition. Therefore, a natural way to qualita- 
tively improve on the above-mentioned constituent quark 
models is to include the pionic degrees of freedom. 

Early investigations of the 7*iVA transition includ- 
ing pionic effects were performed within the chiral bag 
model \100\ HOlj . which was developed as an improve- 
ment to the MIT bag model |102) by introducing an el- 
ementary pion, which couples to quarks in the bag in 
such a way that chiral symmetry is restored 103 . Cal- 
culations within the chiral bag model [104] . found that 
with a bag radius, R, around 0.8 fm one is able to ob- 
tain a reasonably good description for the helicity ampli- 
tudes, as can be seen from the values given in Table [I] 
For such a small bag radius, the pionic effects are crucial 
as they account for around 75 % of the total strength of 
the amplitude A 3 / 2 - The same calculation however yields 
EMR ~ —0.03%, in disagreement with experiment. 

The role of the pion-cloud contributions is also high- 
lighted in Skyrme models [105, 106 , in which the nu- 
cleon appears as a soliton solution of an effective non- 
linear meson field theory. The inclusion of rotational 
corrections in such models, leads to a quadrupole distor- 
tion of the classical soliton solution, yielding a value for 
EMR = —2.3% |106j . consistent with experiment. 

The EMR ratio has also been calculated in models, 
with both quarks and pion degrees of freedom such as 
the chiral quark soliton model (xQSM), which interpo- 
lates between a constituent quark model and the Skyrme 
model |107j . For the two flavor case, one finds EMR = 
—2.1% |108j . fairly close to experiment, considering that 
in the %QSM calculation no parametrization adjustment 
has been made to the N — > A transition. However, the 
magnitudes of the photocouplings, which are given in Ta- 
ble |TJ are largely under-predicted in the xQSM. 

A number of subsequent works have revisited quark 
models, restoring chiral symmetry by including two-body 
exchange currents between the quarks. These exchange 
currents lead to non- vanishing 7*JVA quadrupole am- 
plitudes |109j . even if the quark wave functions have 
no D-state admixture. Such a picture [109 , in which 
the A is excited by flipping the spins of two quarks, 
yields EMR ~ -3.5%, and relates the N -> A and A+ 
quadrupole moments to the neutron charge radius as : 



Q p ^a+ = </ V2, Qa+ = K . (21) 

Using the experimental neutron charge radius, = 
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-0.113(3) fm 2 , Eq. |2p yields : Q P ^ A + = -0.08 fm 2 , 
and Qa+ = —0.113 fm . This value of Q p ^a+ 1S close to 
the empirical determination, given in Eq. (20). In such 
hybrid (quark/pion-cloud) models |110) . the pion cloud 
is fully responsible for the non-zero values of the intrin- 
sic quadrupole moments and hence for the non-spherical 
shape of these particles. As a summary, we list in Table [| 
the jNA photo-couplings A 1 / 2 and A 3 / 2 as well as the 
ratio EMR in various models. 
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TABLE I: Summary of the *yNA photo-couplings A^/2, A3/2, 
and EMR in different models compared with experiment. 



3. Large N c predictions 

Although the results obtained from the different QCD 
inspired models reviewed above may provide us with 
physical insight on the j*NA transition and its relation 
to the nucleon and A shape, they are not a rigorous con- 
sequence of QCD. In the following subsections, we will 
discuss what is known on the 7*A^A transition from ap- 
proaches, which are directly related with QCD in some 
limit, such as the 1/N C expansion of QCD (limit of large 
number of colors) , chiral effective field theory (chiral limit 
of small pion masses or momentum transfers) or lattice 
QCD simulations (continuum limit). 

The 1/N C expansion of QCD [fill Q32] provides an 
expansion with a perturbative parameter at all energy 
scales. This expansion has proved quite useful in de- 
scribing properties of baryons, such as, ground-state and 
excited masses, magnetic moments, and electromagnetic 
decays. For reviews see Refs. [1131 1114] . For example, 
the N — > A transition magnetic moment (J-n^a is re- 



lated to the isovector nucleon magnetic moment as |115j : 
/ i p-i-A+ — (Mp — A*n) /v2 — 3.23 fi ]v , within 10 % of the 
experimental value of Eq. ( |19[). The EMR value was 
shown to be of order 1 /N 2 [116) . Thus its smallness is 
naturally explained in the large N c limit. 

The large N c limit also allows one to relate the A and 
N — » A quadrupole moments via [117] : Qa+ I Qp->A+ = 
2y2/5 + (D(l/N 2 ), Using the phenomenological value 
Qa+ = - (0.048 ± 0.002) fm 2 , which 
1.87±0.08. 

A+ and r 2 was 



of Eq. (|20J yields 
implies Ge2(0) = 

The relation of Eq. (21 1 between Q p . 
also shown |117j to holdin the large N c limit. Further- 
more, it was shown (118] that in the large N c limit : 



EMR: 



CMR = (1/12) R%1 {Ml 



M 2 N )r 2 Jn v , (22) 



with Rna = Mn/Ma, and Ky — k p - K n , the isovec- 
tor nucleon anomalous magnetic moment. Numerically, 
Eq. @ yields EMR = CMR = -2.77%. For EMR this 



prediction is in an excellent agreement with experiment, 
Eq. (20). For CMR, where a direct measurement at the 



real-photon point is not possible, extending the large- N c 
relation to finite Q 2 allows relations with the neutron 
electric FF, which agree well with experiment [1181 1119] , 



4- Lattice QCD and chiral effective field theory 

Lattice QCD provides the possibility of calculating the 
N to A e.m. FFs starting from the underlying theory of 
QCD. The set-up for the lattice calculation of the three 
transition FFs is the same as that used for the extrac- 
tion of the A FFs. The advantage in this case is that the 
connected diagram yields the full contribution. It is cal- 
culated by sequential inversion through the sink for the 
same three simulations as described in the case of the A 
FFs. In addition, recent calculations using Np = 2 + 1 
dynamical DWF generated by the RBC and UKQCD col- 
laborations [120) provide a unitary setup and a further 
check of the results |121j . In all these calculations the 
pion mass is such that the A is still stable. Fig. [TH] shows 
a comparison of the lattice results for G* M1 at the lightest 
pion mass in each type of simulation. There is agreement 
among Wilson fermions and results obtained using the 
hybrid action as well as dynamical DWF. The agreement 
of results using dynamical fermions with the quenched 
results indicate that pion cloud contributions due to pair 
creation are still small at a pion mass of about 330 MeV. 
As also seen for the nucleon e.m. FFs, lattice results un- 
derestimate G* M1 . Chiral dynamics is expected to induce 
large corrections at small Q 2 and such effects can be in- 
vestigated as lattice simulations at smaller pion masses 
become available. 

The CMR and EMR are shown in Fig. [l3j and have 
larger statistical errors due to the fact that G* E2 and G* C2 , 
being sub-dominant, are harder to determine. We show 
quenched results that have the smallest errors as well as 
results obtained in the hybrid action approach and using 
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FIG. 16: Q 2 -dependence of the TV -¥ A FF GJ./i, at the 
lightest pion mass for each type of simulation. Quenched re- 
sults are shown with the asterisks, results with Nf = 2 Wilson 
with the open squares, results using the hybrid action with 
the dotted squares and using DWF with the filled circles. Ex- 
perimental data from Refs. [57J [S3 [70] [73] ED EE2H are shown 
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FIG. 17: The dependence of EMR (upper panel) and 
CMR (lower panel), at Q 2 = 0.1 GeV 2 . The data points are 
from MAMI [68] (blue circle), and BATES [HUE]) (green 
squares). The filled black diamonds are quenched lattice 
results [132] . The open diamond near ~ represents 
their extrapolation assuming linear dependence in m 2 . The 
solid black squares at = 0.36 GeV are hybrid lattice 
results [123] at Q 2 ~ 0.04 GeV 2 (upper points) and at 
Q 2 ~ 0.2 GeV 2 (lower points). The solid red curves are a 
XEFT result [131] . The (blue) error bands represent the esti- 
mate of theoretical uncertainty for the xEFT calculation. 



Np = 2 + 1 DWF QH]. The conclusion that can be 
drawn is that agreement of lattice results on EMR and 
CMR with experiment is better as compared to that of 
G* M1 . Such an agreement is seen also in other ratios, 
indicating that they are less affected by lattice artifacts 
than each of the quantities separately. 

The present lattice QCD calculations are performed 



for quark masses larger than their values in nature. To 
extrapolate to the physical pion mass, one can use the 
XEFT of QCD ptHT25] . xEFT provides a firm theo- 
retical framework at low scales, with the relevant sym- 
metries of QCD built in consistently. The 7*A^A tran- 
sition provides new challenges for ^EFT as it involves 
the interplay of two light mass scales : the pion mass 
and the N — A mass difference. A first study, tak- 
ing into account these two mass scales, was performed 
within the framework of heavy-baryon chiral perturba- 
tion theory [127] . A more comprehensive study was sub- 
sequently carried out [1281 1129j using the "e-expansion" 
scheme. In that scheme, the two light scales in the prob- 
lem: the pion mass e = m 7r /A x sB, with A x ge ~ 1 GeV 
the chiral symmetry breaking scale, and the A-resonance 
excitation energy S = (Ma — Mjv)/A x sb are counted 
as being of the same order, i.e. e ~ S. To allow for 
an energy-dependent power-counting scheme designed to 
take account of the large variation of the A-resonancc 
contributions with energy, the "5-expansion" scheme has 
been introduced |130j . It treats the two light scales e 
and S on a different footing, counting e ~ S 2 , the closest 
integer-power relation between these parameters in the 
real world. It has been applied to the study of the 7* TV A 
FFs (13 lj , and has been used in extrapolating the present 
lattice QCD calculations to the physical pion mass. This 
is shown in Fig. [17] for the EMR and CMR ratios, which 
shows that %EFT predicts strong non-analytic dependen- 
cies on the quark mass for < (Ma — Mjv), invalidat- 
ing simple linear extrapolation in m 2 .. In particular, the 
xEFT results reconcile the lattice results and the rela- 
tively large negative experimental value for CMR. 

For smaller pion masses, where the A becomes an 
open channel, the lattice results will be able to provide 
momentum-dependent phase-shifts [49]. To extract res- 
onance quantities from those will require a fitting pro- 
cedure, e.g. Breit-Wigner or complex pole fits, as done 
when extracting them from experimental multipolcs. 



V. CONCLUSIONS 

In this review, we have presented the experimental 
results and theoretical understanding on the shape of 
hadrons. Although shapes of nuclei have been explored 
over many decades, it is only in recent years that it be- 
came possible to define this question in a theoretically 
rigorous way for hadrons, and perform the experiments to 
answer it. The key concept is to quantify size and shape 
of an extended object through a quantum mechanical 
density operator. For a relativistic bound state system of 
near massless quarks, a probability interpretation is ob- 
tained by considering the system in a light-front frame, 
and projecting its charge density along the line-of-sight. 
We have argued that the resulting transverse charge den- 
sity encodes the information on hadron size and shape. 

On the experimental side, the most accessible and 
best studied reaction to reveal hadron deformation is the 
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N — >• A transition. We have reviewed the state-of-the- 
art experimental techniques, which have allowed to ac- 
curately determine the N — > A quadrupolc amplitudes 
at low momentum transfers, and establish a deformation 
in the N/A system. The quadrupole transitions were 
pinned down on the order of a few % of the dominant 
magnetic dipole transition. A quantitative understand- 
ing of the small, non-zero values of these amplitudes from 
the underlying theory, QCD, is a particular challenge. 
We have provided the historical perspective in which this 
question was addressed from QCD inspired models, high- 
lighting the role the pions play in these transitions. It is 
only very recently, however, that ab initio calculations 
became possible, and state-of-the-art full lattice QCD 
simulations for both the N —> A and A quadrupole FFs 
were able to quantify them. There is an ongoing effort by 
many groups to perform such simulations at pion masses 
approaching the physical value and reducing further lat- 
tice artifacts. 

The theoretical foundations, experimental techniques, 
and lattice QCD simulation methods to access hadron 



deformation through the measurement of quadrupole FFs 
are well established now. We can therefore expect in the 
near future that a refinement of the lattice calculations 
as well as new high-precision experiments with polarized 
beams and polarized targets, will allow to further sharpen 
our understanding of hadron shapes. 
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